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Abstract. The LIGO and Virgo detectors are sensitive to a variety of noise sources, 
such as instrumental artifacts and environmental disturbances. The Stochastic 
Transient Analysis Multi-detector Pipeline (STAMP) has been developed to search for 
long-duration (t>ls) gravitational-wave (GW) signals. This pipeline can also be used 
to identify environmental noise transients. Here we present an algorithm to determine 
when long-duration noise sources couple into the interferometers, as well as identify 
what these noise sources are. We analyze the cross-power between a GW strain channel 
and an environmental sensor, using pattern recognition tools to identify statistically 
significant structure in cross-power time-frequency maps. We identify interferometer 
noise from airplanes, helicopters, thunderstorms and other sources. Examples from 
LIGO's sixth science run, S6, and Virgo's third scientific run, VSR3, are presented. 



Identification of long-duration noise transients in LIGO and Virgo 



2 



1. Introduction 

The Laser Interferometer Gravitational-wave Observatory (LIGO) [1] and Virgo [2] 
detector are part of a network of gravitational-wave (GW) interferometers seeking 
to make the first direct observations of GWs |. These detectors are designed to be 
sensitive to GWs emitted by the coalescence of neutron stars or black holes, core collapse 
supernovae, pulsars, or cosmological processes that happened during the early universe 
[3]. 

In July 2009, LIGO commenced its sixth science run (S6), while Virgo started its 
second and third science runs (VSR2 and VSR3) interrupted between January 2010 and 
July 2010 for a detector upgrade. In this paper, only VSR3 data have been studied. 
LIGO consists of three Michelson interferometers at two different sites. There are 
two located in Hanford, Washington, known as HI and H2, with 4km and 2km arms 
respectively, though only HI acquired data during S6. The other is located in Livingston, 
LA, a 4km interferometer known as LI. Virgo consists of a 3km interferometer located 
in Cascina, Italy, known as VI. 

In order to make confident detections of GWs, it is necessary to understand 
detector noise artifacts. Despite their insulation, the detectors are susceptible to a 
variety of instrumental and environmental noise sources that decrease their detection 
sensitivity [I]. Short, non-astrophysical transient events or glitches can mask or mimic 
real signals. Environmental noise can couple into the interferometer through mechanical 
vibration or because of electromagnetic influence. Seismic motion from human activity 
near the sites, from wind and from ocean waves are among the most common sources of 
these disturbances. Though not as common as ground traffic, airplanes and helicopters 
also affect the GW strain amplitude, or h{t). 

To mitigate the effects of noise artifacts in GW searches, periods of time when 
the detector noise level is unsuitable for GW detection are excluded or vetoed [1]. It 
is necessary to develop techniques to veto non-astrophysical events both effectively and 
safely. A veto is considered safe if no GW signal is dismissed. Vetoes reduce the adverse 
effect of noise on searches for GWs. Probes such as microphones and seismometers 
are placed at sensitive locations around the detectors. Hundreds of channels are thus 
recorded continuously. These are called Physical Environmental Monitoring (PEM) in 
LIGO and auxiliary channels in Virgo. In the following, we will use PEM to refer to 
any of these channels, regardless of detector. A number of veto definition methods have 
been developed using these channels to pinpoint short duration noise transient sources 
[5, 6, 7]. Techniques for identifying narrow spectral features ("lines") that could mimic 
a continuous GW source have been developed as well [S]. 

The Stochastic Transient Analysis Multi-detector Pipeline (STAMP) was recently 
developed to search LIGO and Virgo data for long GW transients whose durations may 
range from 0{s) to weeks. It is a MATLAB implementation of the excess cross-power 

I As presented at the Gravitational-wave Physics and Astronomy Workshop in Milwaukee, Wisconsin, 
January 26-29, 2011. 
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statistic described here [9]. This algorithm uses frequency-time (/t)-maps of GW strain 
cross-power, calculated by cross-correlating strain channels from two spatially separated 
detectors. Pattern recognition methods (including the box- and Radon algorithms) are 
then applied to identify GW signals. A complete discussion can be found in [9]. 

Here we apply STAMP techniques to analyze cross-power /t-maps obtained with 
a PEM channel and h{t) for HI, LI and VI data during S6 and VSR3 runs. By cross- 
correlating h{t) with a PEM channel, there is no sensitivity to astrophysical signals — 
which are not present in the PEM channels — while preserving coherent noise artifacts. 
The same pattern recognition algorithms mentioned above are then used to find noise 
sources. To create the /t-maps, the cross-power is calculated in either 4s or 52s blocks 
over the total time segment requested, depending on the approximate duration of the 
type of noise source. We use overlapping Hann-windowed data. The advantage of 
this method is that only noise sources which couple into the GW strain amplitude are 
revealed. An excess of noise in a PEM channel alone will not induce excess cross-power 
when the PEM channel is cross-correlated with a GW strain channel. By looking for 
structure in these /t-maps (see Fig. 1 for an example), it is possible to identify the 
origin of the noise. 
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Figure 1. /t-map of the cross-power between the GW channel and an accelcrometer 
at the Hanford site using 52s segments. There are strong hnes at multiples of 60Hz, 
corresponding to the harmonics of the 60Hz power mains present at both sites. 

The remainder of the paper is organized as follows. In section 2, we present the 
search for airplane and helicopter noise artifacts. Identification of a noise source covering 
the Crab pulsar frequency during two weeks in S6 is discussed in section 3. In section 
4, we present a search for thunderstorms. Our conclusions are discussed in section 5. 
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2. Airplanes and Helicopters 
2.1. Airplanes 

Airplanes are a source of non-stationary noise for LIGO and Virgo, as they generate 
acoustic noise. An airplane's movement with respect to the detector causes a Doppler- 
modulated signal on an /t-map made with acoustic PEM channels. The airplane's 
frequency (assuming a straight-line, constant-speed trajectory) falls monotonically with 
time. The existing LIGO airplane veto system (called planemon [10]) has been shown to 
flag airplanes observed in microphone channels, and these flags have been shown to agree 
with airplane flight data [9]. A disadvantage to this method is that it does not determine 
if the vibration from the airplanes couple into h{t). To reduce misidentification, it also 
requires multiple microphones to be triggered to identify a signal, and thus planemon 
may miss low amplitude airplane noise that triggers only one microphone. 

In order to use STAMP to identify airplane events, /t-maps of 4s segments 
correlating h{t) with acoustic channels (microphones) are computed in 400s blocks. 
400s was chosen so as to capture the entire airplane track, but not so long that the 
track loses its curvature in an ft-map. The average length for an airplane track during 
S6 and VSR3 was 115s. The frequency range of the /f-map was restricted to 65Hz- 
115Hz, as most of the airplane signal appears in this frequency band. We then take 
the absolute value of the /t-map, as transient noise artifacts can produce complex (not 
positive-definite) cross-power because the cross-power phase depends on the coupling 
of the environmental noise into h{t). Purely horizontal and vertical lines in the maps, 
which can occur due to glitches and instrumental lines, are removed. 

An example airplane signal can be seen in the /t-map on the left of Fig. 2. We next 
apply a Radon transform [11] to each /t-map, an example of which can be seen on the 
right of Fig. 2. The Radon transform produces Radon SNR, which is a weighted sum of 
cross-power along lines parameterized by the impact parameter b, the shortest distance 
to the line from the origin and 6, the angle from the horizontal. Although airplane tracks 
are not quite lines in /t-space, the approximation of the tracks as lines is suitable for 
a simple identification. We then record the maximum Radon SNR, denoted as SNRp, 
which is the value of the largest pixel in each Radon map. 

To determine an approximate threshold for SNRp, a number of /t-maps for 
airplane tracks are checked by eye for each detector site. In Fig. 3, we show a histogram 
of SNRr for the Virgo data. Using histograms like this, we set thresholds on SNRp. 
They are 16.5, 16, 16 for VI, HI and LI respectively. 

For each map with SNRr above the threshold, it is deemed to contain an airplane. 
To fit the track to a line, we set the Radon map to zero except where it is maximal and 
perform an inverse Radon transform back to /t-space (see Fig. 4 for an example). The 
start and stop times of the airplane noise in this frequency band are estimated to be 
the times at which the reconstructed track intersects the edge of the /t-map, and these 
define the veto window size for this airplane. 

This algorithm was run automatically for each day in VSR3 and S6d, which occurs 
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Figure 2. Left: /i-map of the cross-power between the GW channel and an acoustic 
channel. The airplane track can be seen on the left side of the plot. Right: Radon 
map of this event, where the impact parameter b is the shortest distance to the line 
from the origin and is the angle from the horizontal. The bright spot on the right 
panel corresponds to the airplane's signal in Radon space. 
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Figure 3. Histogram of SNRr for above threshold airplane events (in red, dotted 
outline) , near threshold events (green dashed outline) and bclow-thrcshold background 
events (in blue, solid outline) in a week of Virgo data. The dotted black line corresponds 
to the threshold of 16.5. The red, green and blue distributions are arbitrarily 
normalized for the purpose of plotting. About 93% of the maps correspond to non- 
airplane noise, 6% arc airplanes, and 1% arc ambiguous. Entries with SNRr ~ 25 
record the maps with SNR > 25. 
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Figure 4. Left: Radon reconstruction of the airplane event in Fig. 2. Right: the 
"veto" window created with the reconstruction plot in Fig. 2. 



Detector 


Time Analyzed 


Number of Analyzed Microphones 


Number of Airplanes per Day 


Dead Time 


Vetoed HW Injections 


HI 


S6d 


13 


28 


0.08% 


4/9752 


LI 


S6d 


9 


21 


0.14% 


9/8464 


VI 


VSR3 


4 


2.9 


0.09%. 


0/3083 



Table 1. Results from the airplane analysis for VSR3 and S6d. For each detector, 
the time period analyzed, number of analyzed microphones, number of airplanes per 
day, dead time, and vetoed HW injections are listed. Numbers of airplanes per day 
detected have been cross-checked with other monitors. 



in the last two months of S6. Each day was divided into 400s blocks, and h{t) at each 
detector was correlated with the available microphones. For each block deemed to 
contain an airplane, the GPS times of the identified airplane were saved. For each day, 
a list of the times of the identified airplanes was compiled. 

Over the entire run, to estimate the effect of airplanes in the detectors, the 
deadtime (percentage of data excluded by the airplane flag) was calculated. In order to 
determine if the veto is "safe," we use simulated GW waveforms of burst and coalescing 
binary inspiral signals injected into the interferometer by applying a force on one of the 
mirrors of the Fabry- Perot cavity in one of the interferometer's arms to make sure that 
the veto is not issued by a GW event [G]. These results are listed in Table 1. For all three 
detectors, the number of vetoed injections is consistent with statistical expectations for 
a safe veto. Even so, a number of airplanes were found to be accidentally coincident 
with hardware injections. Due to the way vetoes are currently used, these periods would 
have been vetoed. But STAMP is able to determine the frequency band of the noise 
disturbance which can be well separated from a real GW signal frequency. It is possible 
to observe a GW event in the presence of an airplane (as long as the GW and airplane 
events do not overlap in the /t-map). It therefore seems plausible that during future 
science runs, this capability will be implemented to veto less data from GW analysis. 
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2.2. Helicopters 

During VSR3, helicopters were a noise source for the Virgo detector, with tracks similar 
to those of airplanes, although their tracks contained more visible harmonics (see Fig. 5 
for an example). Some of the flybys were close enough to Virgo to degrade detector 
sensitivity dramatically. Due to their similar noise structure, the airplane algorithm 
identified a number of helicopter events which were vetoed. 



I 
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Figure 5. Left: Spectrogram of a microphone in tire central building of the Virgo 
detector at a time when a helicopter was flying near to the interferometer. The arrows 
indicate instances of helicopter noise. Right: /t-map of the cross-power between the 
GW channel and this microphone at this time. This plot shows that the noise of the 
helicopter is coupling with h(t). 



3. Identification of line near the Crab Frequency 

The Crab pulsar, which has a GW signal frequency of approximately 59.55Hz, is an 
important target for LIGO and Virgo's continuous GW searches [12]. During two 
weeks of LlGO's S6 run, an initially unidentified noise source appeared in h{t) near 
the Crab frequency, /t-maps correlating h{t) and all available PEM channels, using 
52s segments, were calculated for this time. The noise was observed in a number of 
seismometers located inside of the central building, where the beam splitter and two of 
the test masses reside. A chilled water booster pump in the mechanical room was found 
to be vibrating at this frequency, coupling into h{t). A number of weeks before, the 
pump's load had been reduced, causing its frequency to approach the Crab frequency. 
Subsequently, the load was returned to normal, eliminating the noise, as can be seen in 
Fig. 6. 
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Figure 6. Left: /i-map of the cross-power between the GW channel and an 
accelerometer near to the water pump on June 8, 2010 (before the normal load was 
resumed) . Please note that the line was not stationary and so the noise does not appear 
continuously in the /t-map. Right: ft-m&p of the same on June 13, 2010 (after the 
normal load was resumed). 

4. Thunderstorms 

Thunderstorms are a source of magnetic and acoustic environmental noise. They can 
prevent the detectors from running stably. In order to identify thunderstorms, we create 
/t-maps correlating the GW channel and microphones, as well as maps of the GW 
channel and magnetometers. We then search for broadband events of high cross-power 
in these maps. The times of these events are recorded for both the microphones and 
magnetometers, and events coincident in both the microphones and magnetometers 
within a 10s window are flagged as thunderstorm candidates (see Fig. 7 for an example). 
During S6d, 7 events were identified at Hanford, while 26 were identified at Livingston. 
During VSR3, 1 event was identified at Virgo. 

To test the efficiency of the algorithm, the identified events were compared with 
206 identified lightning events one night at HI using a camera on the roof. During 
this time, only 2 of these lightning hashes were identified by the algorithm, resulting 
in a low efficiency of 1%. To ensure that these were not accidental coincidences, we 
gave the coincident events artificial time offsets, or "time slides," looking for lightning 
events. Time slides between -lOOs-lOOs in 20s intervals were performed, resulting in 
zero accidental coincidences. This suggests that the algorithm is detecting actual 
lightning flashes, but at a very low rate. There are a few possible reasons for this. 
One is that it is unlikely that each lightning strike has enough strength to appear in a 
microphone, magnetometer, and the interferometer strain amplitude. Also, it is likely 
that instead of detecting the magnetic flelds of the lightning, we are detecting magnetic 
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fields from power system glitclies due to ligliting strikes tliat liit power lines or other 
power distribution facilities. This was confirmed by observing events in voltage monitors 
on site coincident with the lightning events. 




time (s) time (s) 

Figure 7. Left: /t-map of the cross-power between the GW channel and a microphone 
during a coincident event at HI. Right: /t-map of the cross-power between the GW 
channel and a magnetometer at that time. 

5. Conclusions 

We have realized many noise hunting applications with the Stochastic Transient Analysis 
Multi-detector Pipeline (STAMP), and others are still in development. STAMP data 
quality flags were generated for use by LIGO and Virgo inspiral [13, 11] and burst 
[15, 16] GW searches. Looking ahead towards advanced detectors, STAMP will be an 
important tool for identifying not only the noise sources currently affecting the detector 
but also the new noise sources that will inevitably appear. 
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